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^ ■ We examine effects of the next-nearest-neighbor repulsion on electronic states of a one- 

dimensional interacting electron system which consists of quarter-filled band and interactions 
of on-site and nearest-neighbor repulsion. We derive the effective Hamiltonian for the elec- 
trons around wave number i/cp {k-p: Fermi wave number) and apply the renormalization group 
method to the bosonized Hamiltonian. It is shown that the next-nearest-neighbor repulsion 
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O ' makes 4fcF-charge ordering unstable and suppresses the spin fluctuation. Further the excitation 

gaps and spin susceptibility are also evaluated. 

^ ' KEYWORDS: organic compounds, charge density waves, spin density waves, quarter-filling, long-range Coulomb 

00 ■ . , ,. . 

^— N ' interaction, spin gap 

0\ 
O 

o 
o 



§1. Introduction 



't^ . Quasi-one-dimensional organic conductors, many of which are the 2:1 charge transfer salts, exhibit 

\^ • a variety of physical properties depending on temperature, pressure and so on.El^ A one-dimensional 

'^ ' (1"D) interacting electron system at quarter- filhng is a basic model for understanding such exotic 

Q , properties. The electronic states of these materials have been often theoretically investigated by 

use of the model with only the on-site interaction [/ (> 0). Such a model with the on-site repulsion 
1^. has been used to explain spin density wave (SDW) states observed experimentally. 

Vh ' Recent experiments on these materials indicate a coexistence of SDW with charge density wave 

(CDW).ii&i In (TMTSF)2PF6, 2A:f-SDW state coexists with 2kF- and Akp-CDW below the SDW 
transition temperature, a' while the two kinds of CDW's disappear below 3-4 KB' The material 
(DI-DCNQI)2Ag shows the magnetic order of 2A;f-SDW below 5.5 KE> However, the analysis of 
^^C-NMR spectrao suggests the 4A;f-CDW ordering below 220 K, which is much higher than the 
transition temperature of the magnetic order. Such a 4A;f-CDW order has been also observed in 

* E-mail : yoshioka@phys.nara-wu.ac.jp 



X-ray studyB) It is known that the coexistence of SDW and CDW observed experimentally cannot 
be explained by a model with only the on-site repulsion. Therefore, it is necessary to treat models 
with the long-range components of Coulomb interaction, which may take crucial roles for these 
coexistent states. 

Effects of several long-range interactions have been theoretically studied especially for a 1-D 
quarter-filled interacting electron system, and the rich phase diagram has been found. Numerical 
diagonalization of the model with both U and the nearest-neighbor interaction, Vi,Q'EfoE3' shows 
that the insulating phase appears for large strengths of both U and Vi , and that the superconducting 
state becomes the most dominant fluctuation for large Vi and small U . The mean-field theories have 
predicted the some coexistent states of SDW and CDW.tirllSllalliP For the system with U and Vi, a 
transition occurs from a pure 2A;f-SDW state to a coexistent state of 2A;f-SDW and 4A;f-CDW with 
increasing Vi.EJ^ It is maintained that such a coexistent state gives rise to the charge ordering in 
(DI-DCNQI)2Ag.H' The transition has been examined by evaluating the commensurability energy 
corresponding to the S/sp-Umklapp scattering.H' Kobayashi et al. found that the next-nearest- 
neighbor repulsion, V2, results in a coexistence of 2kF-ST)W and purely electronic 2kF-CDW x3> 
The problem has been reexamined by Tomio and one of the present authorsEJ) and it has been 
found that the coexistent state of 2A;f-SDW and 2A;f-CDW is followed by 4A;f-SDW but without 
4A;f-CDW. The next-nearest-neighbor repulsion is the most promising candidate for the origin of 
the coexistence observed in (TMTSF)2PF6.i'i 

In the previous works, we have analytically derived the 8/cF-Umklapp scattering for the electron 
system at quarter-filling with U and Vi, and investigated the electronic properties by applying 
renormalization group (RG) method to the bosonized model.E3'ES* The phase diagram has been 
determined on the plane of U and Vi, which is qualitatively the same as that derived by the 
numerical approach. We found that the insulating state for both large U and Vi exhibits the order 
of Akp-CDW with the most dominant fluctuation given by 2A;f-SDW. Such a state is similar to the 
coexistent state predicted by the mean-field theory and seems to correspond to the experimental 
observation in (DI-DCNQI)2Ag that the characteristic temperature of the charge ordering is much 
larger than that of 2A;f-SDW. 

In the present paper, we extend the previous methodtjUB) to the system with U, V\ and V2, and 
clarify the role of V2 for the electronic states and the excitations.llZP We find that V2 makes ^k^- 
CDW state unstable and suppresses the spin fiuctuation. The phase diagram is determined and 
compared with that by the mean-field theory. Electronic states are further examined by evaluating 
the excitation gaps of the charge and spin fluctuation, and spin susceptibility. 

The plan of the paper is as follows. In §2, the effective Hamiltonian is derived and expressed in 
terms of bosonization with the phase variables. The phase diagram and excitations are determined 
by utilizing RG method in §3. Based on the RG equations, the spin susceptibility is calculated in 



§4. Section 5 is devoted to summary. 

§2. Model and Formulation 

2.1 Effective Hamiltonian 

We consider a 1-D electron system at quarter-filling. The Hamiltonian is given hy 7i = TCo + Ti-mt, 
where 

= J2i^K - ^i)a^K,a'^K,a, (1) 

Ka 

Wint = X] 1 ^'^CT'-o-"'i,fTnj>' + VlUj^arij + l^a' + V2nj>nj+2,a' \ 
jaa' ^ ^ > 

X SKi+K2-K3-K4,GaKi,aO'K2,(T''^K3,a'aK4,cr- (2) 



Here t and fi denote a transfer energy between the nearest-neighbor sites and a chemical potential, 
respectively, ex = — 2tcosi^a with a lattice constant a, and — vr/a < K < vr/o. The quantity 
aj^(= 1/\/Nl J2k G-'dp^ 0-) i^ ^ creation operator of the electron at the j'-th site with spin 
(j(= lb), Uj^^ = a- fjCLj^a-, G = and ib27r/a, and Nl is the number of the lattice. 

In order to obtain the effective Hamiltonian for the states near zt/cp (fcp = 7r/(4a)), the one- 
particle states are divided as shown in Fig. 1, i.e., dk~,a = o,K,a for —vr/a < K < — 7r/(2a), 
Ck,-,a = CiK,cr for — 7r/(2a) < K < 0, Cfc,+,o- = aK,cr for < i^ < 7r/(2a), and dfc.+.o- = ^^",0- for 
7r/(2a) < K < n/a, where A: is deviation of the wave number from ztfep for Cfc^±^o-, and itS/cF for 
dk,±,a- Hereafter the one-particle states expressed by dk^±^a- (cfc.i.o-) are called as the upper (lower) 
band. In terms of 0^,^,0- and d^^p^a (p = i), W is written asTi = Tio + Y^i=Q Ti-int,!, where 

T-CO = ^[{(^pkp+k - M)Cfc,p,^Cfc,p,a + (eSpfep+fc - fJ')dk,p,ad'k,p,aJ , (3) 

pkcr 

Wint.O = jYJ2Yl y^^'''-'' ~ ^20) Ck+q,p,a(^k'-q-p,a''^k\p,a'Ck-p,a 
k,k',q V u,a' ^ ^ 

A:,fc',g P cr.cr' ^ ^ 

+ 7 IZ mZ ( ^'^'^'.-'^ + ^1" + ^20) 4+g,p,a4'-g,p,<x'Cfc',p,<T'Cfe,p,a, (4) 

Wi„t,l = 7 IZ im (^«'^'t',-'^ - 2^20) 4+g,p,a4'-g,-p,a'Cfc',p,<7'Cfe,p,(r 
k,k',q V a,a' 

+ tJ2J2J2 {Ua5„,^-„ - 2V2a) 4+g,p,a4'-g,p,a'Cfc',-p,a'4,p,a, (5) 

k,k',q P a,a' 




^K 



Fig. 1. The energy dispersion, ek = —2t cos Ka, in the present system. 



k,k',q P a,a' 

k,k',q P a,a' 

+ 7 H im (^Sa',-a - V2a] {4+g,p,a4'-g,-p,a'4',-p,a'4,p,a + h.C.} 
k,k',q P a,a' ^ ' 

+ tJ2 J2J2(^^^^'^'^ + ^^^^ + ^^2^)4+q,P,adk'~q,~pydk',-p,a'Ck,p,a 
k,k',q P a,a' 

+ jJ2J2J2 [^Sa',-a-Via + V2aj {cl^q^p^A'-q,-p,a'dk',p,a'dk,-p,a + h.C.} 
k,k',q P a,a' ^ ^ 

+ 7 H im (f^a<5a',-a - 2Fia + 2V2a) 4+g,p,<x4'-g,-p,a'Cfc',p,a'4,-p,a 

k,k',q P a,a' 

+ 7 H im ( -Y^'^''-'' ~ ^1" + ^^a) {ci+g,p,A'-g,p,a'dk',~pydk,-p,a + /i.e.} 

A:,fc',g P cr.cr' ^ ^ 

+ 7 II im (^a(5^',-a - 2lia + 2V2a) 4+g,p,<x4'-g,p,a'Cfc',-p,a'4,-p,a 

+ 7 II mZ (^«<^'T',-'T - 2^2a) cl^q,p,a4'-q,p,a'dk',-p,a'Ck,-p,a, (6) 

fc.fc'.q P cr.cr' 

Wint,3 = 7 IZ mZ (^"'^'^'.-'^ - 2^2a) A+q,p,adl'-q,-p,a'dk',p,a'dk,p,a 



k,k',q P a,a' 



+ 7 IZ Y.Y.(^^^^''-^ - 2^2a) 4+g,p,a4'-g,p,a'4',-p,'^'Cfc,P,a, 



L 



(7) 



k,k',q P a,a' 



k,k',q V cr,cr' ^ ^ 



+ J mini 'Y^^''-'' + ^^" + ^2") 4+g,p,a4'-g,-p,a'4',-p,a'4,p,a 



In the above expressions, L = Nlu and the matrix elements of the interaction are calculated at 
K = ib/cF for the lower band and at i^ = ib3/cF for the upper band. The elements, {Ua/25at -„ + 
Via + V2a), {Ua/25cri -a — V2a) and {Ua/25a> ^^a — Via + V2a) are originated from the interaction 
processes with the momentum transfer, 0, 7r/(2a) and vr/a, respectively. Since there is no Umklapp 
scattering for only the lower band, i.e., the first term of Hq a-nd Hjnt.O) it is crucial to take account 
of the effects of the upper band in order to describe the insulating state of the quarter-filled 1-D 
electron systems. 

By averaging the contribution from the upper band, the effective Hamiltonian describ- 
ing the states near it/cF is derived as follows. Representing Ck^p^a and d^^p^a in terms of 
Grassmann variables, the partition function corresponding to eqs. (3)-(8) is written by Z = 
/ ^[4,p,a4,p,a]^[4,p,^Cfc,p,^]e-'5, where 

4 

5* = S'oc + S'od + 2^ 5*1111,1, (9) 

f/3 



fP 

"^Oc = 51 / dr4 p ,^(a^ - /i + epfep+fe)cfc,p,a, (10) 

SQd=z_\ / drd^ (3r - /i + e3pfcp+fe)c^A:,p,(T, (H) 

^ Jo 



p,k,a ' 
Smt,i = / drHint.j, (12) 

Jo 
with f3 = 1/T (T: temperature). After integrating with respect to dk,p,a, the partition function, 
Z, is given as Z = Zod/P[4p,^Cfc,p,a]e"^<=«, where Zq^ = 1 1^[dl^p,adk,p,a]e~^°'', and S^s = Soc - 
ln{^exp(-X;f=o'S'int,i)) ((-4)^ = Z^^J jV[dlp^^dk,p,a]e^^"'^A). The perturbative expansion of ^cfr 
leads to the effective Hamiltonian. 

Up to the second order, the normal processes, Hint,!!, which denote backward scattering and 
forward scattering, are derived as 



~q,~p,cr'^k',p,a'Ck,-p,a 



Wint,n = jJ2J2J2\ (9l±S^~a' + 9l\\Saa')cl+g,p^„cl,^ 
k,k',q V u,u' I 

+ (fi'2±'^a-a' + 5'2||'^a<7')4+(?,p,<T4'-g,-p,a'Cfe' -p,a'Cfc,p,a 

+ (ff4±<^a-a' + 54||'^a<7')4+g,p,<74'-g,p,a'Cfc',p,a'Cfc,p,a \- (13) 



9.. = '^-V,a-ADJ'^-V,a\('^-V,a + V,a]^ 



The coupling constants are given as follows, 

5-111 = -V2a - ADi{-V2a){-Via + V2a 

Ua 
= --y.a-4Z^,^2 -VV2 

(72|| = Via + V2a - 2Di {-V2af - 2Di {-Via + V2af , 
52± = ^ + Via + V2a - 2Di (^ - V2a] - 2Di (^ - Via + V2a^ 
5411 = Via + V2a - 2D2 {-Via + V2af , 
<?4± = ^ + Via + Vsa - 2D2 l-^ - Via + ^20 



(14) 
(15) 
(16) 
(17) 
(18) 
(19) 



where Di = {S-Kta)-^ V2ln{V2 + l) ~ 1.25/(87rta) and D2 = {87rta)-^2V2/{V2 + l) ~ 1.17/(87rto). 
Note that eqs. (14)-(17) satisfy the condition, g2± — 5211 + ^iH = 9i±j which corresponds to SU{2) 
symmetry. The diagrams which lead to eqs. (14)- (19) are shown in Fig. 2, and the detailed cal- 
culation is given in Appendix A. The 8A;F-Umklapp scattering is obtained from the third order 
expansion shown in Fig. 3. The Hamiltonian for the Umklapp scattering is written as 
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Fig. 2. Diagrams for the normal scattering. Here the dashed and wavy hnes express the Green function of the upper 
bands and the interaction, respectively. 



Wi/4 = 91/4^2 h^XXp,+i'-p,+f{Xp~i^-p-f, 



(20a) 



51/4 = ^{{Ua-2V2a)\Ua-Wia + V2a)+i{V2af{Ua-2Via + 2V2a)}, (20b) 
where V'p.o- = ^/yLJ^k^ '^k,p,c7- Detailed derivation of eqs. (20a) and (20b) is given in Appendix 
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Fig. 3. Diagrams for the Sfep-Umklapp scattering. The dashed and wavy hnes express the Green function of the 
upper band and the interaction, respectively. 



B. Thus, the effective Hamiltonian is written as 

HcS = Woc + Wint.n + W1/4, (21) 

where Hqc = Y^pkai^vkp+k - M)4,p,aCfc,p,a- 

2.2 Bosonization 

We Hnearize the energy dispersion as epkp+k — H ~ pvpk with vp = \f2ta and utihze the bosoniza- 
tion method. The phase variables for the charge (spin) fluctuation, Qp and (/)p, {Qf, and ^a) are 



defined as 

vri 



q^O ^ k,p 



(22) 



9^0 ^ fe,p 

where a~^ is the ultraviolet cutoff. Those fields satisfy the commutation relations, [6u{x), <piy'{x')] = 
i7rsgn(x — x') 61^1,1 {u, u' = p 01 a). The electron operator is expressed as Vp.o- = (27ra)^^/^ 
exp{i{p/2){ep+p(j)p + aea + crp(l)a)} exp(i7rHp^) with Hp+ = p/2J2p' ]^p'+ and Hp_ = p/2J2p' Np'^ + 
J2p' ^p'+ where Np^ = / dxV^lo-V'po- • In terms of these phase variables, the effective Hamiltonian is 
expressed as Weflf = "^p + "^o- , where 

np = ^Jdx l^idjpf + Kpid^cPpA + ^^iZi^ |dxcos40p, (24) 

n^ = ^Jdxi^^{dJ^)^ + K^{d^<t>^f^ + ^^Jdxcos29^. (25) 



Hereii'^ = y/B^/A^ and v^ = v^^/ByAy with Ap(^) = 1 + {54|| +(-)94±+fi'2|| + (-)52±-5i||}/(™f) 
and Sp(^) = 1 + {5411 + (-)54± - 52|| - (+)52± + fi'i||}/(7rt'F)- The coupling constant, 51/4, is given 
by 5i/4 = 5'i/4/(2'?i"a)^- In terms of the phase variables, the order parameters for the 2/i;f-SDW, 
2fcF-CDW, 4A:f-CDW and AkY-hond order wave (4A:f-B0W) are written as follows, 

02fcp-SDW = E ^e-'^P'^F^V^^^.p,, ex sin(2fcFX + Op) sin 6^, (26) 



per 



02fep-cDW = E « ''^^^^'''^iai'-P,^ « cos(2A:fx + 0p) COS e„, (27) 

per 

04fcp-CDW = E e-'^^''^"V'J,+V'J,-^-p,-V'-P,+ « cos(4A:fx + 26 p), (28) 

p 

04fcp-B0W = E e"''^'^^"(ip)V'J,+V'J,-V'-p,-V'-p,+ oc sin(4A;FX + 26 p), (29) 

p 

where the order parameter of 4A:f-B0W is given by the 4A;F-component of nj+nj+i^-. Note that 
the third order perturbation for ln(exp(— X]i=o 5'int,j))d shown in Fig. 3 generates the other 8A;f- 
Umklapp scattering term, cos 40p cos 20o- ) which mixes the charge and spin fluctuation. Since the 
scaling dimension of this term is smaller than that of cos A9p or of cos 26 ^j , we discard the term. 
The effects of the term are discussed in §5. 

§3. Excitations and Phase Diagram 

We investigate the excitations and the electronic states in the limit of low energy by using the 
RG method. At first, we consider the spin degree of freedom. By introducing / = ln(a'/a) as the 
new length scale a' (> a), the RG equations for Ti^r are given as follows ,113' 

^^K.{l) = -]^Gl^{l)Kl{l), (30) 

j^Gi^{l) = [2-2K^{l)]Gi^{l), (31) 



written as^^" 



where the initial conditions are given by K„{0) = K^j and Gi_l(0) = 2gi^/TTV(j. The behaviors 
of Kfj{l) and Gi_l(/) are as follows. When gi± > (i.e., U/2 > V2), the quantities, Gi_l(/) and 
Kfj{l) tend to and 1, respectively, and the excitation becomes gapless. Then the long-range 
correlation functions are given as (sin0(j(x) sin^(j(0)) ~ x^-'^ln ' {x) and {cos 6 (j{x) cos 6 cr{^)) ~ 
x~^ ln^^'^(x).ll3' For U/2 < V2, the spin gap appears due to Gi_l(/) -^ —00 and K(j{l) — > 0. In 
this case, the phase, ^o-, is locked as mod vr at the low energy limit and then {svii9a{x)sm.9(j{{))) 
vanishes exponentially and {cos9fj{x) cos9fj(0)) remains finite at the long distance. Therefore 2kp- 
SDW does not exist. Thus it is found that the spin fluctuation is suppressed by V2. 
Next we investigate the charge degree of freedom. The RG equations corresponding to Tip are 



±Kp{l) = -8GJ/,(l)K'^{l), (32) 

|Gi/4(/) = [2-8Kp(/)]Gi/4(/), (33) 

with the initial conditions given by Kp{0) = Kp and ^1/4(0) = giu/{2iTVp). Correlation 
functions are calculated from the solutions of eqs. (32) and (33) as, {sin 6 p{x) sin 6 p{0)) = 
{cos9p{x)cos9p{0)) oc exp{- Jq''''^'^ dlKp{l)}, {cos29p{x)cos29p{0)) oc exp{- /J'"'^/" d/[4Kp(/) + 
2Gi/4(/)]} and (sin20p(x) sin20p(O)) oc exp{- /d""'/" d^[4Kp(/) - 2Gi/4(/)]}.0) Equations (32) and 
(33) are calculated analytically asEj) 

Gl/^{1) - l/{2Kp{l)) -2\nKp{l) = A, (34) 

where ^4 is a numerical constant determined by the initial conditions. Equation (34) is obtained 
by rewriting eq. (32) as G\i^{l) = {l/S)dK-^{l)/dl and Kp{l)Gli^{l) = -{l/8)dlnKp{l)/dl, which 
are substituted into eq. (33). By choosing some parameters, we show eq. (34) on the plane of 
Kp and G1/4 in Fig. 4. In the region (a), the non-linear Giu term becomes irrelevant leading to 
metallic state. Here the correlation functions with long distance are given by {sin9p{x) sin9p{0)) = 
{cos9p{x)cos9p{0)) ~ x-^p^°°'> and (cos26lp(x) cos26'p(0)) ~ (sin 26lp(x) sin 26*^(0)) ~ x~^^''^°°\ 
In particular, the initial values located between the thick solid curve and the thick dotted one are 
renormalized to 1/4 < Kp{oo) < 1/3 and ^1/4(00) = 0. Both regions (b) and (c) corresponds to the 
strong coupling regime with insulating states, in which |Gi/4(/)| — > 00 and Kp{l) -^ 0. In the region 
(b) where Gi/4{1) — > —00, the correlation functions, (cos20p(2;) cos20p(O)), (cos0p(x) cos6'p(0)), 
and (sin0p(x) sin0p(O)) remain finite at the long distance, whereas (sin26'p(x) sin20p(O)) decays 
exponentially. Then both 2A;F-density waves and 4A;f-CDW states are stabilized. On the other 
hand, in the region (c) where Giu{l) -^ 00, one finds the finite value of (sin20p(x) sin20p(O)), 
(cos0p(x) cos6'p(0)) and (sin0p(x) sin6'p(0)) but the exponential decay for (cos26'p(x) cos20p(O)). 
Therefore 4A;f-B0W and 2A;F-density waves are stabilized and 4A;f-CDW is completely suppressed. 
In Fig. 4, the closed circles, open circles and closed squares denote the initial values for Vi/t = 1, 




Fig. 4. Scaling flows of eq. (34) on the plane of Kp and Gi/4JiHtZP The thick solid curves denote the boundaries 
for metallic regions (a) and insulating regions (b) and (c). The fixed point for the thick dotted curve is given by 
{Kp,Gi/4) = (1/3,0). The closed circles, open circles and closed squares denote the initial values of Vi/t = 1, 3 
and 5, respectively, with the fixed U /t = 6, where respective points from right to left axe V2/t — 0, 0.6, 1.2, 1.8, 
2.4 and 3.0. 



3 and 5, respectively, with the fixed U/t = 6, where respective points from right to left are V2/t = 
0, 0.6, 1.2, 1.8, 2.4 and 3.0. Here a ~ 2a/7r is used-llB* Since increase of V2 with fixed U and Vi 
moves {Kp, G1/4) to the region (c), it is found that the next-nearest-neighbor repulsion V2 stabilizes 
4A;f-B0W states. 

By combining above consideration on the charge and spin degrees of freedom, we obtain the 
phase diagram, which is shown on the plane of Vi/t and l^/t with fixed U/t = 6 in Fig. 5. Here 
the most dominant state of the respective region is as follows; 2/i:F-SDW metal (I), 4A;f-CDW metal 
(Ha), 4A;f-B0W metal (Hb), 4A;f-CDW insulator (HI), 4A:f-B0W insulator (IV), and 2A;f-CDW 
-|- 4/cf-BOW insulator (V). The region (I) with small Vi/t and 1^/t, which exhibits the gapless 
excitation for both charge and spin, denotes the metallic state with 2A;f-SDW as the most dominant 
fluctuation. In the region (HI), the insulating ordered state of 4A;f-CDW appears, where 2A;f-SDW 
exists as the most dominant fluctuation. This comes from the fact that the parameters located in 
the metallic region (a) in Fig. 4 move to the insulating region (b) with increase of Vi and fixed small 
V2, whereas the spin excitation remains gapless. Note that between the insulating state of 4:kp- 
CDW and the metallic state of 2/cf-SDW, the metallic state with the most dominant fluctuation of 
4A;f-CDW exists (region (Ha)). Such a region is obtained for 1/4 < Kp{oo) < 1/3 and gi/4 < (see 
Fig. 4). With increasing V2/t, the parameters move from the region (b) to the region (c) in Fig. 4. 
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Fig. 5. The phase diagram on the plane of Vi/t and V2/t at U/t = 6 where the thin dashed curve expresses gi^4 — 0. 
The states of respective regions are given by 2fcp-SDW metal (I), 4fcF-CDW metal (Ila), 4fcF-BOW metal (lib), 
4feF-CDW insulator (III), 4fcF-BOW insulator (IV), and 2fcF-CDW + 4A:f-BOW insulator (V). In the region (III) 
and (IV), 2A;f-SDW exists as the dominant fluctuation. 



In this case, the 4A;f-CDW state disappears and the ordered state of 4A:f-B0W appears (region 
(IV)), where the 2A;f-SDW state remains as the most dominant fluctuation as long as the spin 
excitation is gapless. The metallic state of 2/i;F-SDW in the region (I) turns into the insulating 4:kp- 
BOW state (region (IV)) when V2 increases with fixed small Vi as seen from the closed circles with 
Vi/t = 1 in Fig. 4. With increasing V2 further, the spin gap appears for U < 2V2 and the insulating 
2A;f-CDW + 4A;f-B0W state is realized in the region (V) since the correlation functions of 2kp- 
SDW and of 4/i:F-CDW are exponentially decayed at the long distance. In contrast to the region 
(Ila), there exists the metallic 4A:f-B0W state (region (lib)) corresponding to 1/4 < Kp{oo) < 1/3 
and Qiu > between the region (I) and (IV). 

Figure 6 shows the charge gap Ap (a) and the spin gap A^- (b) as a function of V2 with choices of 
Vi/t = 6, Vi/t = 5 and Vi/t = 3.5 for U/t = 6. By using the solution of the RG equations, (30)- 
(33), the quantities Ap and Ao- are determined by the following condition, Ap = VpO"^ exp(— /p) 
with \Gi/^{lp)\ = l)13 and Ao- = iio-a"^ exp(— /g-) with Gi±{lfj) = —2. Note that the spin gap 
determined by such a method shows quantitatively good agreement with that obtained from the 
exact solution^) for Ao-/t < 1 in case of negative U Hubbard model. For Vi/t = 6 and 5, the charge 
gap is finite except for the boundary between the region (III) and the region (IV) in Fig. 5. In 
the regions (IV) and (V), Ap increases monotonously with V2. In the region (III), Ap for Vi/t = 5 
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Vp/t 



Fig. 6. The charge gap Ap (a) and the spin gap A^ (b) as a function of V2 at U/t — 6 with choices of Vi/t = 6 
(sohd curve), Vi/t = 5 (dotted curve) and Vi/t = 3.5 (dashed curve). 



does not show a monotonous behavior as a function of V2, whereas Ap for Vi/t = 6 decreases with 
increase of V2. Such a behavior of Ap for Vi/t = 5 may come from the fact that the boundary 
between the region (Ila) and the region (III) is not a monotonous function of V2. As mentioned 
above, the condition for appearance of the spin gap is given hy U < IV2 and then independent of 
V\. However, as seen in Fig. 6(b), the magnitude of the spin gap does depend on the value of Vi, 
which enhances the spin gap. 

Here we note the case of small U . The metallic state is realized for all regions corresponding to 
Fig. 5, i.e., the region for V\ < U and V2 < U . In this case, the new metallic state appears for 
U < 2V2 in which the spin gap exists as shown in just below eq. (31). In this metallic state, the 
most dominant fluctuation is 2A;f-CDW for Kp{oo) < 1. 

Finally, we consider correspondence between the present state and the coexistent state found 
by the mean-field theory. As discussed in the previous paperJlS the insulating Ak-p-CDW states 
in the region (HI) corresponds to the coexistent state of 4A;f-CDW and 2A;F-SDWciP because the 
most dominant fluctuation of 2/i:F-SDW exists in the region (HI). In the insulating 4A;f-B0W state 
(IV), the fluctuations of both 2/cf-CDW and 2A;f-SDW develop, whereas the 4A:f-CDW state is 
completely suppressed. Therefore, the state of the region (IV) may correspond to the coexistent 
states of two kinds of 2A;F-density wavesll3 without 4/cf-CDW.ES 



12 



§4. Spin Susceptibility 

In this section, the spin susceptibihty of the present system is investigated. We utihze the formula 



developed in ref. 21, which is effective at relatively low temperatures since it is based on the RG 



analysis using the linearized band. From eqs. (30) and (31), one obtains 



5i±(0 



9i± 



l + 25i±(vr7;0)-ir 



(35) 



where gi±/{TTVa) ^ 1 and v'^ = vp <1 + {g^ — g4±){TTv-F) ^ >■ By using eq. (35), the spin suscepti- 
bility, Xcr(T), is calculated asElP 



Xa{T) 

xliT) 



xm 



2l-{g4±-g4\\+gi±mx°p{T)' 

2 \-^ ,. fi^pkp+k) — fi^pkp+k+q) 



- — > lim 



^ ^ 9^0 epkp+k — ^pkp+k+q 



(36) 
(37) 



where /(e) = l/(e('-'')/'^ + 1) and / = \n{vYa-^/T). In Fig. 7, x<7(r)/x°(0) is shown where x°(0) 
denotes the susceptibility in the absence of the interaction at T = 0. At high temperatures, 
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Fig. 7. Normalized spin susceptibility, Xo-(r)/x'^(0), where x'^(O) is the susceptibility in the absence of the interac- 
tion at r = 0, as a function of T/t for the several choice of V2 with fixed U/t = 6 and Vi/t — 5. The black and 
white circles denote the values of T = for U > 2V2 and U — 21^2 , respectively. Inset : Spin susceptibility at T = 
as a function of V2/t with fixed U/t = 6 and Vi/t — 5. 



the spin susceptibility increases with decreasing temperature. In the limit of zero temperature, 
the susceptibility approaches to a constant with infinite slope for U > 2V2 as seen in Hubbard 
modehElP Since it is due to logarithmic temperature dependence of 5'ij_(/), such an anomaly is not 
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observed for U = 2V2- In case of the strong coupling, i.e., U < 2V21 the susceptibihty vanishes 
near the temperature, Tsg = 'Vfo:~^ e^p{TTv'^/{2gi±)}, which originates from the fact that we 
use the solution of the weak coupling expansion, eq. (35). It is expected that the susceptibility 
for T ^ Tsg behaves as [Aa/T)^'"^ ex.p{—Aa/T}c3' where Ao- is the spin gap. The inset shows 
Xo-(0)/x'^(0) as a function of V2/t with fixed U/t = 6 and Vi/t = 5. The interaction V2 suppresses 
the spin fluctuation monotonously even for the region without the spin gap. On the other hand, 
Vi enhances X(t(0)/x°(0). These are due to the effect from the upper band because Xa{0)/x^{0) = 
{1 + (5411 - 5'4±)(vrvF)~^} with 5411 - g^^ = -Ua/2{1 - 2D-2{Ua/2 - 2Via + 2V-2a)}. Note that 
the spin susceptibility in case of Vi = V2 = is different from that derived in ref. ^ because our 
formula includes the virtual process to the upper band. Spin susceptibilities both with and without 
the virtual processes are compared with that derived from the exact solutioncfP in Appendix C. 

§5. Summary 

We investigated the electronic states and excitations of the one-dimensional electron system 
at quarter-filling with on-site (C/), nearest-neighbor (Vi) and next-nearest-neighbor (V2) repulsive 
interactions. By taking account of the effects of the upper band in a perturbative way, the effective 
Hamiltonian has been obtained in terms of the phase variables and analyzed by the renormalization 
group method. 

The following roles of V2 have been found. The spin fluctuation is suppressed with increasing 
V2 and the spin gap is obtained for V2 > U /2. Such a characteristic behavior can be seen in the 
spin susceptibility, X(t{T)-, where Xo-(O) decreases with increasing V2 and vanishes for V2 > U/2. 
Note that the decrease of Xcr(O) by V2 for U/2 > V2 is due to the effects of the upper band. In 
addition, the interaction suppresses the 4A;f-CDW state and stabilizes the 4A;f-B0W state. As 
a result, the insulating state of 4/cf-CDW realized for large U and Vi changes into the insulating 
state of 4A;f-B0W with increasing 1^. The insulating 4/i;F-B0W state moves to the insulating states 
of 4kp-JiOW + 2A;f-CDW with increasing V2 further. Thus the charge distribution changes from 
(O o 0°) of 4kF-CDW to (O O °°) of 2/cF-CDW with increase of V2 where Q (°) denotes the rich 
(poor) charge at the respective site. The excitation gaps were determined as a function of V2 with 
several choices of Vi at U/t = 6 from the solutions of the RG equations. In the states of insulating 
4A;f-B0W and 4A;f-B0W + 2A;f-CDW, the charge gap increases as a function of V2, whereas the 
gap is not a monotonous function of V2 in the insulating 4A;f-CDW state when parameters are 
located near the boundary of metal insulator transition. The spin gap appears independently of 
the value of Vi. However, the magnitude of the spin gap does depend on Vi, which enhances the 
spin gap. For the case of weak on-site repulsion, the ground state for Vi < U is metallic. When V2 
is added to the state, the metallic state with spin gap appears, where the dominant fluctuation is 
given by the 2A;f-CDW. 

In the present treatment, we have neglected the other S/cF-Umklapp scattering term. 
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COS 40p COS 2^0-, which mixes the charge and spin degrees of freedom and has the scahng dimen- 
sion less than that of cos20o- or of cos4^p. The importance of such a term has been pointed out 
for half-fihing extended Hubbard model with U + VifBE3c3' in which a term cos 2^p cos 20o- exists 
besides cos 26 p and cos 2^0- (the definitions of the phases and the interaction correspond to ours). 
In the half-filling model with U + Vi, the coefficient of cos 29 p and of cos2^o- vanishes at U = 2V\, 
whereas the term cos26pCos26(^ is finite. Thus the new phase appears near U ~ 2Vi.EZP In the 
present model, the coefficient of cos 46 p and that of cos 26„ become zero simultaneously at the 
point, U = Vi = 2V2- In this point, the coefficient of cos 4^p cos 20o- also vanishes (see Appendix 
B). Therefore, it is expected that such a term does not change the phase diagram qualitatively, but 
may result in only a quantitative change. 

The appearance of the spin gap and the variation of the charge order with increasing V2 are of 
interest and may be observed by the increase of the pressure along the conducting chain.B3) 
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Appendix A: Derivation of the normal processes 

From Fig. 2, the matrix elements are given as 

51,0-0-' = -'^2fcF,o-o-' — 4— 2^ 2^X2fcp,o-o'-'^4fcF,o-o-'G'p(ie„, /c)G_p(— ien, — /c), (A-1) 

'^ -7r/(4a)<fc<7r/4a '^n 

92,aa' = Xo^aa' " 2— ^ ^ {^2%,oo' + ^Ifcp.aa' } ^^(ien, k)G-p{-ien, -k), (A-2) 

' -7r/(4a)<fc<7r/4a ^n 

94,00' = ^0,00' - 2— Yl YxIkF,aa'G-p(.ien,k)G^p{-ien,-k), (A-3) 

^ -7r/(4a)<fe<7r/4a ^n 

where gi^^a = gi\\ and gi,a-a = gi± (i =1,2,4). The quantities, Xo,<jo', X2kp,aa' and X4kp^aa' are the 
matrix elements of the interaction, Xo^o-o-' = Ua/25cr-a' + ^lO + ^2^, X2kY,au' = f^o/2(5o-,-o' — ^2^ 
and X4fcp o-o-/ = Ua/25fj^^fj' —Via + V2a, respectively. The green function of the upper band is given 
by Gp(ie„, k) = l/(ie„ — e^pkp+k + /^)- The particle-particle propagaters are calculated as, 

^1 = «^ H ^Gp(ie„,fc)G_p(-ie„,-/c) 

-7r/(4a)<fc<7r/4a <5n 

= 7 J2 -Ku -v{l-2/(e3p%+fc-/^)} 

dy- ^^ = -^ln(V2 + l), (A-4) 



S-Kta Jo smy + l/\/2 Sirta 
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-^2 = -^ Yl ^Gp(ie„,fc)Gp(-ie„,-fc) 

'"■'^ -7r/(4a)<fc<7r/4a ^n 

1 1 

1 /-Va 1 1 2^2 ,, ,, 

dy- ^ —^ = -— ^— , (A-5) 



4:TTta Jo sin y + cos y + \/2 Sirta ^/2 + 1 

where /(e) = l/{e^'^ + l). Here f {(-^pkY±k — lA is treated as zero because the energy scale we consider 
is sniaUer than —^ = v2t. Equations (A-l)-(A-5) give rise to eqs. (14)-(19). 

Appendix B: Derivation of Umklapp scattering 

As shown in Fig. 3, the S/cp-Uniklapp scattering consists of the interaction processes, in which the 
right moving four electrons are scattered into the left moving states and vice versa. Such processes 
are derived from (•S'j^^^^ 6*1111,2) /2 among the third order perturbation. The processes in Fig. 3(a), 
(b) and (c) are written as Sa, Sb and Sc, respectively, where 

1 fP rl3 rl3 

^'^ = Ti^HYl J2 / '^'^l'^^2kp,aia[ Yl Yl / dr22X2fcp,^2^^ ^ ^ / dnX^^^^^.^^, 
P ai,(T[ ki,k[,gi 0-2,0-2 «2,fe2 '92 03,0^ fe3,A;?5,(j3 

'^ \ (4i+Qi,-p,oi(n)4;-Qi,-p,oi(^i)cfe;,p,o;(n)4i,-p,oi(n) 

^ 42+g2-p,02('^2)Cfc^_g2_p,^^(T2)Cfc/^^p^^/(r2)42,-p,02(^2) 

'^ 43+93 -P,03(^3)4^_g3,-p,^^(r3)cfc/^,p,4(r3)cfe3,p,^3(r3))d + h.c. L (B-1) 

1 /•/3 rl3 rf3 

^b=JjYY Y J ^Tl'^^2kp,a,a[ Y Y J ^^'i'^X^k^ ,a2CT'^ Y Y J ^'^l'^^ ik^ ,^Z<t'., 
P ai,a[ki,k[,qi 02,02 A;2,fc2:92 03,03 fc3,fc3, 53 

X (4i+gi,p,oi(n)4j_<;,_p,^^(n)cfej,p,ai(n)cfei,p,ai(ri) 

^ 42+g2-p,02('^2)Cfc^_g2_p,^^(r2)Cfc/_^p^^/(r2)42,-p,02(^2) 

X 43+93 -p,03(^3)4;(-Q3,-p,o^(^3)Cfc^,p,<7;,(T-3)43,p,03(T3))d, (B-2) 

sc=^YY Y jyn2x,,^^,^,,Y Y jyr22x,,^^^^,,Y Y jyn2x,,^^^,^^, 

P o-i,a[ki,k[,qi 02,02^2,^:2,92 03,03 fc3,fc3,ij3 

X (4i+gi,p,oi(n)4;-gi,-p,o;(^l)Cfe;,p,o;(Tl)Cfei,p,ai(ri) 

X 42+g2-p,o2('^2)4^-g2,-p,o^(^2)Cfc^,p,4(r2)42,-p,02(^2) 

^ 43+93 -P,03 (^3)4^-93, _p,a,;(^3)4f,,p,o;,(^3)Cfe3,p,a3(T3))d- (B-3) 

In eq. (B-1), the average, (• • •)^ is calculated as, 

- \4i,-p,oi (ti)42,-P,02 (^2)43+93 -P,03 ('^3)4^,-93, -p,o;, ('^3))^ 

= -G_p(ti - Ts, ki)G-p{T2 - T3, ^2) ('^fc3+g3,fc2'^^-g3,fci<^03,02<^o^,oi - ^k3+q3,ki^ki^-q3,k2^(T3,(TiSai^^a2) 

- -^'^('^1 - T3)6{t2 - T3) [Sk3+q3,k2^k'^-q3,ki^'^3,'r2Sa'^,ai " 43+<?3,fcl 4;j-<?3,fc2 '^03,01 (^4,02) • (B-4) 
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Here the green function, Gp{T,k) = j3 ^Y^e^^ '''"'' Gp (ie„, fc), is replaced by -l/(2^/2t)(5(r). By 
substituting eq. (B-4) into eq. (B-1), 



8t2 

f/3 



P cri,a[ 0-2,0-2 



1 fP 

^ JJ ^ j drAp^,(gi -g3)^p^2(<?2+g3)^poi(-9i)^po^(-92), (B-5) 

where ylp^o-(g) = J2k^k+q p a'^k,-p,a- In a similar calculation, one obtains S^ = Sa and 

"^c = -^ ZJ zJ Z2 '^^2kp,aia[ X 2X2fcp^^j^^ X 2X2kp,a3ai 
P iTi,a[ 02,0-3 

1 fl^ 

^ 73 51 / dr^P'^,(gi-g3)ylp^/(^2 + g3)^po^(-9i)^po3(-^2)- (B-6) 

91,92,93-^° 

Here we introduce the field operator, ■(/'p,o = ^/\^J2k^''^^'^k,p,a and write the Hamiltonian corre- 
sponding to the above action, Ti-'^u, as 

P cri,(j[ CT2,02 
+ o72 2J zJ zJ 2X2fcp^o^o^ X 2X2A;p,<jio(, X 2X2A;p,030i / dxAp^^Ap^'^Ap^'^Apa-^, (B-7) 

p 01,0^02,03 

where Ap^o- = '4^ta'^-p,u- The above Hamiltonian can be rewritten as 'H'^u = W1/4 + H^^ + HJm, 
where 

Wl/4 = TT2 Z^ I (2-^2fcF ,0-0) (2X4fcF ,00 + 2X4fcp 0— o) + (2X2fcF ,00) 2X4fcp CT- o| / ^^^pa^p-cr 

+ o72 2J^^2A:f,oo(2X2A:f,o-o) j<^xAp^Ap_^, (B-8) 

"^1/4 ~ 77? / ,{2(2X2fcp^o— o) (2X2fcF, 00) (2X4fcF, 00 + 2X4fcp o— o)} / <^xAp^Ap^cr 
p,o -^ 

+ ^ E {6^2fcF,o-o(2X2fep,,,)2 + {2X2kp,.-.f] J dxAl^Ap^^, (B-9) 

^5/4 = ^ E {(2^2fcF,oo)'(2X4fep,..) + ^(2X2fcp,..)3| I dxA^^. (B-IO) 

In the bosonized form, eqs. (B-8), (B-9) and (B-10) are proportional to /dxcos40p, 
/ dx cos 49 p cos 2^0- and / dx cos 49 p cos 4^0-1 respectively, and the scaling dimension of the respective 
term is given by 2 — 8Kp, 2 — 8Kp — 2K^ and 2 — 8Kp — 8K(^. Therefore eq. (B-8) has the largest scal- 
ing dimension. By substituting 2X2fcp ^0-0' = U aS^r^^a' — 2V2a aiid2X/ik^,j„i = Ua6f^^^a-' — 2Via+2V2a 
into eq. (B-8), one obtains eqs. (20a) and (20b). We note that the coupling constant of the Hamil- 
tonian (B-9) is proportional to 2X2fcp, cr-o = Ua — 2V2a. Equation (B-10) becomes important for 
i^CT < 1/4 in the insulating state. However, the term, cos 2^0- is already relevant for K^^ < 1, i.e., 
U < 2V2- Then we can safely neglect Tt 



1/4- 
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Fig. 8. Spin susceptibility of Hubbard model at T = as a function of U/t, where X°(0) is the susceptibility in the 
absence of the interaction. The solid (dotted) curve expresses the spin susceptibility with (without) the virtual 
process to the upper band shown in Fig. 2. The solid circle is derived from the exact solution in ref. E4|. 



Appendix C: Spin Susceptibility at T = for Hubbard Model with [/ > 

In case of Vi = V2 = 0, C/ > and T = 0, the spin susceptibility, eqs. (34) and (35), normalized 
by x'^(O) reduces to, 

Xa(0) 1 



X°(0) 



1 



Ua 

2'KVp 



[I - UD2) 



(C-l) 



We note that the L'2-terni originates from the virtual process to the upper band, which is not taken 



account in ref. 21. In Fig. 8, we show Xcr(0)/x''(0) with and without the virtual processes by solid 
curve and dotted curve, respectively, which are compared with the exact one ( closed circle ).ESAs 
seen in Fig. 8, the spin susceptibility with the virtual processes is very close to the exact solution 
compared to that without the processes. It turns out that the virtual processes are necessary for 
quantitative estimation of the spin degree of freedom. 
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